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Abstract 



A relativistic two-body wave equation, local in configuration space, is 
derived from the Bethe-Salpeter equation for two scalar particles bound by 
a scalar Coulomb interaction. The two-body bound-state wave equation is 
solved analytically, giving a two-body Bohr-Sommerfeld formula whose en- 
ergies agree with the Bethe-Salpeter equation to order a 4 for all quantum 
states. From the Bohr-Sommerfeld formula, along with the expectation val- 
ues of two remaining small corrections, the energy levels of the scalar Coulomb 
Bethe-Salpeter equation are worked out to order a 6 for all states. 
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I. INTRODUCTION 



For two-body atomic systems such as the hydrogen atom and positronium, it would be 
a great advantage to have the energy levels to order a 4 given by a relativistic two-body 
wave equation local in configuration space and soluble analytically. We would like to have 
a two-body counterpart to the one-body Coulomb Dirac equation |IJ with its exact solution 
as found by Gordon and by Darwin in 1928. 

This is a very old problem. In this paper we solve it for a system that is analogous to 
real atoms but simpler: two scalar particles of masses m and M, bound by a scalar Coulomb 
potential —a jr. 

The approach is to start from the Bethe-Salpeter equation Q for the system. A Bethe- 
Salpeter equation is used because in QED it is the standard two-body bound-state equation, 
known to be true, and in which higher-order corrections are well understood. In QED the 
binding interaction in the Coulomb gauge kernel is —oc/r. 

In the present paper we will treat the scalar Coulomb system's Bethe-Salpeter equation 
as the true equation. We will require that the energy levels of any reduction of the Bethe- 
Salpeter equation agree with the levels of the original Bethe-Salpeter equation. 

We reduce the Bethe-Salpeter equation to a two-body relativistic bound-state wave equa- 
tion which is local in configuration space. The locality of this reduction is in contrast to the 
Salpeter reduction Q which has been the standard for the last half century. The Salpeter 
reduction contains operators such as \J— V 2 + m 2 which are non-local in configuration space 
and difficult to calculate with. 

The relativistic two-body Coulomb wave equation has an analytic solution. The energy 
levels are expressed by a two-body Bohr-Sommerfeld formula. 

The Bohr-Sommerfeld formula predicts the correct energy levels to order a 4 . They agree 
to order a 4 with the Bethe-Salpeter energy levels. 

The formalism also allows calculation of the Bethe-Salpeter energy levels to order a 5 
and a 6 . Because the levels up to first relativistic order, a 4 , are given by an exact solution, 
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only first-order perturbation theory is needed for the level calculations up to order a 6 . 
This is in contrast to the standard Salpeter reduction, in which the first relativistic order 
(a 4 ) already needs first-order perturbation theory, while higher orders would need second- 
order perturbation theory. In the present paper, explicit expressions for the Bethe-Salpeter 
equation's energy levels up to order ct 6 are given for every quantum state. 

In Section II of this paper the Bethe-Salpeter equation for the scalar system is written 
down, its Salpeter reduction is reviewed, and its energy levels to first relativistic order (a 4 ) 
are calculated from standard perturbation theory for later comparison. 

In Section III the local relativistic two-body wave equation is derived from the Bethe- 
Salpeter equation. It is proven that the wave equation's first-order relativistic corrections 
are exactly the same as those of the originating Bethe-Salpeter equation for any static scalar 
interaction, so that to calculate relativistic corrections in first-order perturbation theory it 
is just as accurate to use the relativistic two-body wave equation as the Salpeter reduction. 

In Section IV, the relativistic two-body Coulomb wave equation is solved analytically and 
a two-body Bohr-Sommerfeld formula for the energy levels is derived which does predict the 
energy levels of the Bethe-Salpeter equation correctly to order a 4 . 

In Section V the formalism is used to calculate the Bethe-Salpeter energy levels up 
to order a 6 . One set of corrections is given by a simple expansion of the Bohr-Sommerfeld 
formula to order a 6 . The other two correction terms are calculated in first-order perturbation 
theory. 

To our knowledge the work here is the first analytic solution of a relativistic reduction 
of a Coulomb Bethe-Salpeter equation accurate to order a 4 . It may also contain the first 
accurate two-body Bohr-Sommerfeld formula. We have also demonstrated that corrections 
to the energy levels of a Coulomb Bethe-Salpeter equation up to order a 6 can be calculated 
in first-order perturbation theory, once the solution is found analytically to order a 4 . 

It is hoped that techniques like these may one day lead to similar results for real atomic 
systems such as positronium and the hydrogen atom. 
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II. SCALAR BETHE-SALPETER EQUATION AND SALPETER REDUCTION 



A. Bethe-Salpeter Equation 

The notation will be the following. Let the masses of the bound particles be m,M. Let 
the mass of the bound state be denoted by E, while defining the bound-state wave number 
below threshold as (5, as well as the particles' individual CM bound-state energies as t and 
T, in the following way: 

E = ym 2 - f3 2 + y/M 2 - p 2 (1) 

2 i A™2 Ayf2> 



t,^^= E ' + %- M ' , (2) 

r^ ^ 11 ^ (3) 

The CM energy-momenta of the particles are written (p, t + p ), (— p, T — p ). 
The Bethe-Salpeter equation for the bound-state vertex function T will be 

r = I * S T (4) 

with 

p 2 + m 2 -2e-(t+p ) 2 p 2 + M 2 -ze- (T-p ) 2 ^ 
The scalar interaction kernel is written 4mM/(k 2 ), or 4mMI(r) in configuration space. The 
star denotes an integration over 4-momentum with a factor \ji present. The factors 2m,2M 
have been moved into the numerators of the propagators of the scalar constituent particles 
in order to make the dimensions of the scalar functions the same as in spin-| systems. 

B. Salpeter Reduction 

Defining 0(p) = J dp°S(p,p )T(p,p°)/27ii, equation (^) when multiplied by S and inte- 
grated over p° gives the non-local Salpeter equation equivalent to (f|): 



p 2 + m 2 + Jp 2 + M 2 + Z(p)I(r) 



[r) = E<j>{v) 



in which Z(p) is a correction operator on the interaction J(r): 



2mM 



Vp 2 + m 2 + VP 2 + M 2 



Vp 2 + mVp 2 + M 2 VP 2 + m 2 + VP 2 + M 2 + E 



(6) 



(7) 



C. First-order Relativistic Correction to Energy Levels 



To calculate the first-order relativistic corrections to the Schrodinger energy levels pre- 
dicted by equation @ the non-relativistic Schrodinger equation is as usual defined to be 

~2 



m + M + — + I(r) 
2fi 



fr (r) = ^o^o (r) 



(8) 



in which fi is the reduced mass. The non-relativistic bound-state wave number fi will be 
defined by 



E = m + M - (3 2 /2/2 



(9) 



Then the Salpeter reduction (H) is expanded in powers of p 2 and (3q in the standard way, 
giving the correction AE = E — Eq\ 



AE 



l n (3 2 i /V fi 2 " 

8 \r^ + W^) + 2^ mim + YA^ 2 + W 2 



(10) 



This constitutes the first-order relativistic correction for the energy level of the original 
Bethe-Salpeter equation (|J). 
In the scalar Coulomb case 



Ur) = - 



a 



(11) 



equation ( |T0D shows that the energy levels of the scalar Coulomb Bethe-Salpeter equation 
have this correction to order a 4 : 

,2 \ 



AE. 



(4) 



a 4 /i ( jj? ji 2 \ a 4 fi I ^ /r 



y Coul iV3(2L + 1) \m 2 M 2 j 8N A \ mM J 
where L is the angular momentum and N is the Bohr quantum number. 
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III. RELATIVISTIC WAVE EQUATION 



A. Derivation 

We return to the original Bethe-Salpeter equation @ and reduce it to a three-dimensional 
form in a manner different from that of Salpeter. 

Because S peaks sharply at p° = when p 2 is not too large, it can be approximated 
Hby 

4mM 2mS(p°) 
g ° (p ' p) = ~ 2E J^+W) ( 3) 

The error of the approximation is defined by 

R = S - S (14) 

The Blankenbecler-Sugar [KJ correction interaction U is defined by the equation 

U = I + I*RU (15) 

Then from the Bethe-Salpeter equation (§) it is easy to deduce the equation T = U * S T. 

Next, defining 0(p) by So T = 27u5(p°)0(p), it follows that, with the notation U = U(p; q), 

,o o . . . 4mM r da . ,, . , . 

(p 2 + /3 2 )0(p) = -— J j^y 3 U(p, 0; q, O)0(q) (16) 

Equation (|i~6|) is exact, in terms of the original Bethe-Salpeter equation (f|). 

For a static kernel /(k 2 ), equation (TE) for Z7 simplifies, because in its solution by iteration 
the integration over the relative-energy variable only acts on R. After that integration, 
equation (|T5| ) becomes 



U(p, 0; q, 0) = /[(p - q) 2 ] + J ^/[(p " 1) 2 ]^(1 2 )^(1, 0; q, 0) (17) 

in which 



/l9N mM 1 1 . 

^ 1 = -PT S , r + , r } 18 

£ | VI 2 + m 2 (VI 2 + m 2 + t) VI 2 + M 2 ( VI 2 + M 2 + T)^ ' 

where the energy-dependent constants t and T were defined in (||) and (|3]). 



B. Local Equation 



Equation (|TT| ) for U leads immediately to a local wave equation accurate to first rela- 
tivistic order. 

Because in many systems the momentum is small compared to the constituent masses, 
it is a reasonable starting point to approximate (|TB|) by its value at l 2 = 0: 



R E = R(Q) 



mM 
~E~ 



1 



+ 



f 



(19) 



m(m + t) M(M + T) 

Then in configuration space equation ( JT7| ) can be iterated once to give a local approximation 
to U: 



U E (r) = I(r) + R E I 2 (r) 



(20) 



Defining in addition the quantity 



m + M 
E 



equation (|T6|) gives 



{p 2 + 2fiZ E [l(r) + R E I 2 (r)] } 0(r) = -/5 2 0(r 



(21) 



(22) 



which is a two-body wave equation local in configuration space. Equation ( |2"2"D contains 
energy-dependent constants, instead of the non-local operators of the Salpeter reduction 
(|). The eigenvalue (3 2 determines the energy E through equation ([[]). 



C. Equivalence Proof 

Equation ( p2[) can be expanded around the Schrodinger equation (§) to calculate the 
energy corrections given by (p2|) to the first relativistic order. 

Expression (|l|) for the energy must be expanded to fourth order in /3 , so that 

AE = "I + w) ~ 2> < 23 > 

7 



In evaluating A/3 2 = /3 2 — /?q from fl22|), it is only necessary to express Re to zero order in 
j3 , as (/i 2 /m 2 + /i 2 /M 2 )/2/i. In addition Z E need only be expanded to second order in /3 , 
as 1 + /3g/2mM; in the I 2 term Ze can be taken to be 1. 

From these values (^) immediately gives the first-order correction 

a/? 2 = — rXf( J )-f4 + ^)( /2 ) ( 24 ) 

m + M \m z M z J 

Substitution of (|24| ) into (|23| ) shows that the expression for AE" is the same as the expression 
(p!0|) derived by means of the Salpeter reduction. 

This constitutes a proof that the local relativistic two-body wave equation (|22|) correctly 
gives the first-order relativistic energy corrections to the Bethe-Salpeter equation (Q). A 
similar proof for spin-| particles was derived some time ago 0. 

D. Infinite-Mass Limit 



In the limit M — > oo, equation (p2|) becomes 



jp 2 + 2m/(r) + -^J 2 (r)| 0(r) = - (m 2 - t 2 ) 0(r) (25) 

where t, defined in (H), is the bound-state energy of the lighter particle. 

To first relativistic order the term 2ml 2 /(m + t) can be replaced by I 2 . Then ( ^5|) 
becomes 

{p 2 + [m + /(r)] 2 }0(r) =t 2 0(r) (26) 
which is the customary Klein-Gordon equation for a scalar particle in a fixed field. 

IV. SCALAR COULOMB SOLUTION 

Specialising to the scalar Coulomb interaction (pT|), the Blankenbecler-Sugar correction 
potential given by ( |TTD becomes to second order 

^Goul(0 = -"+^ (27) 









jp 2 + 2^Z E 




) 







and the relativistic two-body wave equation fl22|) becomes 

= -0 2 <f> (28) 

The radial reduction of equation ( p8|) has the same singularities as that of the Coulomb 
Schrodinger equation, and can be solved the same way. For angular momentum L the radial 
component can be expanded for radial quantum numbers n = 0, 1, 2, • • • as 



e-*G9r) e 5>G9r) 

i=0 



(29) 



which as usual gives a recurrence relation between the coefficients dj. The existence of do, 
and the termination of the series, require the conditions 



< \i I L + \) +a 2 2nZ E R E -~ 



(30) 



and 



n + l + e 



(31) 



Then from ( |5T| ) and expression (|I]) for the energy E, it is easy to deduce the Bohr- 
Sommerfeld formula 



E 



\ 



m 2 + M 2 + 2mM 



\ 



n + l + e 



(32) 



Expanding in the usual way, with Z E and R E inside e only needed to zero order, it is found 
that the Bohr-Sommerfeld expression ( |32|) does predict the correct a 4 correction (|T2| ) to the 
energy levels of the scalar Coulomb Bethe-Salpeter equation, as expected. 



V. ENERGY LEVELS TO ORDER a 6 
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A. Introduction 



Up to this point it has only been shown that the reduction of the Bethe-Salpeter equation 
to a relativistic two-body wave equation as described in the previous section reproduces the 
results of the standard Salpeter reduction to order a 4 , albeit through an analytic solution 
instead of perturbation theory. One way of investigating whether the present reduction may 
be more useful than the method of Salpeter would be to find out whether higher-order terms 
in the energy levels can be calculated more easily than the Salpeter reduction allows. In 
the present section we will calculate the energy levels of the Bethe-Salpeter equation 
analytically to order a 6 for all quantum states, using only first-order perturbation theory. 

With the definition 

Ai?(l 2 ) = R{\ 2 ) - R E (33) 
the expansion of equation (|T7|) for U becomes (with / standing for —a/r) 

U = U Coul + IARI + IRIRI + IRIRIRI + ■■■ (34) 
We will evaluate the contributions to the energy levels up to order a 6 due to: 

• f^Qoul' wnose levels are contained in the Bohr-Sommerfeld formula ([^); 

• I ARI, which will have an a 5 contribution for L = 0, as well as a 6 contributions; 

• the order-a 6 contribution of the rest of the series, 

U 6 = IRIRI + IRIRIRI + ■■■ (35) 
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B. Bohr-Sommerfeld Formula 



The Bohr-Sommerfeld expression (|32|) can be expanded further to get the sixth-order con- 
tribution to the energy of the Blankenbecler-Sugar potential Uq qu ]_, equation (p7|). Defining 
the dimensionless constant 



fi 2 /i 2 



m 



the sixth-order energy contribution is 



^^Bohr-Somm a ^ j 64iV 6 



12c - 5c 2 - 1 3c 2 c 2 1 

+ 8iV 5 (2L + 1) ~ 2iV 4 (2L + l) 2 ~ iV 3 (2L + l) 3 J ^ ' 

We recall that L is the angular momentum and N is the Bohr quantum number. To obtain 
this expression the constants Zg, -Re and e are expanded to the order in a needed. 



C. Correction IARI 



Following the steps outlined in Subsection A above, we calculate the a 5 and a 6 energy 
corrections due to the term I ARI in (33), where AR is defined in (p3|). The correction is 



(6) 



IARI 



(38) 



Up to order a 6 it is sufficient to use non-relativistic wavef unctions. Use of the Schrodinger 
equation (H) then gives 

1 [dp 



AP (6) 

AR ~ (2^ J (27T) 3 ' 



'olPJ 



p 2 + /3 2 A J R(p 2 )0 o (p) 



(39) 



To evaluate the integral fl39|) for every quantum state, it is convenient to change 
coordinates from p-space to the surface of Schwinger's unit sphere [0] in Cutkosky's 4-space 
PJ. The transformation is (with i = 1, 2, 3) 



6 



p 2 + /? 2 



P 2 -/3 2 

P 2 + /3 2 



(40) 
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(Henceforth (3 is written for (3q.) The polar coordinates on the unit sphere are conventionally 
denoted by (0i, 9 2 ,4>) M where 2 , 4> are the usual angles of 3-space, and 



x = cos Q\ = £4 



(41) 



The element of surface area is 



dQ= (sin0i) 2 (sin0 2 )d0id0 2 £ 



(42) 



and dfl and dp are related by 0] 



dtt 



2(3 



dp 



(43) 



_P 2 + (3 2 _ 

Surface harmonics on Schwinger's unit sphere are denoted by Y NLM , where N — 1 > 
L > \M\ > 0. The quantum numbers L,M have their usual meaning in 3-space, and 
iV = 1,2,3,... is the Bohr quantum number. As usual they are related byiV = n + L + l 
where n — 0, 1, 2, . . . is the radial quantum number in p-space. 

The standard representation of the normalised surface harmonics is [Ef 



>W(0i,02,0) = 2 L+1 r(L + i 



n\ N 



(l-x 2 )tC n L+1 (x)F LM (0 2 ,0) (44) 



27rr(A^ + L + 1) 

in which Ylm(02, <P) is normalised on the surface of the unit sphere in 3 dimensions, and the 
coefficient is determined by the orthonormality relation || 

nT(N + L + l) 



(1 - x 2 ) L+ ^C^\x)C^\x) dx = S m , n 



(45) 



n! N \T(L + 1)] 2 2L+1 

The momentum-space eigenfunctions 0o(p) are proportional to (1 — x) 2 Ynlm- With 
the normalisation requirement 

dp 



! 



(27T)* 



>olPj 



(46) 



we find that 



!>olPj 



2tt 
J 



[1 — x) 2 Y NLM (0i, 2 , ■ 



(47) 
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Then substitution of equation ( f47|) into (|39|) gives 



{6) _ 2^n\N[T{L + l)f ^ 



vrr(A^ + L + l) /x 2 ' 



(48) 



in which 



+ 1 



dx (1 — x ) 



2 A#(p 2 ) 
1 — X 



(49) 



Although Ai?(0) = 0, it is not possible to expand Ai?(p 2 ) to first order as a Taylor series 
in p 2 , because the integrals in (|3~9"D and (|49|) would diverge when L = 0. In order to evaluate 
for all states including L = it is necessary to express A_R(p 2 ) exactly. 
With successive definitions, referring to equations (|40|) and ([41]): 



2/J 2 
t 2 



.4 



2^ 

y2 



(50) 



n , , x / Vy + V^T^ t 2 (m + 2t) I 



(51) 



F(x) = 



it can be shown that 



mM f 1 



E 



{^0(1 - x, a, t, m) + ^D{1 - x, A, T, M) 



(52) 



Vl — x cfe 

Expression (|49D can now be integrated by parts. To the required order, it becomes 

87V 



(53) 



<L,n 



—a S 



L.O 



3/i 



/i 3 /i 3 



(54) 



in which 



J, 



L,n 



[ +1 dxF(x)4~ {VTTxil - x 2 ) L 
J-i dx I 



c^Hx) 



(55) 



Note that the 5l,o term is of order a, not a 2 . It will give a contribution of order a 5 fi to the 
energy. In atomic physics the Coulomb potential gives an a 5 contribution when L = 0. In 
the present model and formalism that term appears in the AR correction. 
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The remaining integral ([55]) stays convergent as a, A — > inside the square roots in the 
functions D contained in F(x). To the required order that limit may be taken, and fl55|) 
becomes to lowest order 

1 



Jl,u = -ot 



/iiV 2 



4 4 
jjL fl 



(56) 



in which 



K 



K 



(i) 

L,n 
(2) 



-1 d 

dxVl ~ X— < \/l + x(l - x 2U 
4 J-i dx 



L+l, 



i 2 



X 



+1 W 1 

ax—== 
i vT= 



d 



X 



dx 



VT+~x(l — a; 



2\L 



C: 



L+l, 



(57) 
(58) 



The expression for if^i can De integrated by parts for all L. The expression for K^' n 
can be integrated directly by parts for L > 1, while for L = the integration by parts 



may be done by subtracting C„ +1 (l) 



from 



<^ +1 (x) 



in the integrand beforehand. The 



subtracted part is calculated separately and added back after the integral is evaluated. With 
that understanding, we have 



K 



(i) 

L,n 



3 

+- 



-i a 



X 2 ) L+ 2 



C^\x) 



K 



(2) 
L,n 



3 f+ 1 0_ 
'4 7-1 



1 



-dx 



x 



X 



-dx 



x 



(59) 
(60) 



To evaluate these i"Ct, )n integrals it has been necessary to derive the following equation 
H (in which m < n): 



dx- — 

-l 



x 



2\LA 



tC™{x)C™(x) 



z — X 



0Fr(iV + L + l) 



+1 (1 _ f2^n+L+ 



dt- , 
i (z - t) n+1 



2 JV+L r(L + i)r(A^ + |) 

With m = n the limit z — *■ 1 gives 

_ 3 v^r(L + |)r(iv + L + i) 



(61) 



(62) 



8T(n+ 1)T(L + l)T(2L + 2) 
To evaluate ( |60"D one can differentiate (|6T| ) with respect to z and take the limit z — > 1. 
For L = the subtracted integral is calculated before taking the limit 2 — > 1; the value is 
then finite. For all L, including L = 0, the result is 
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K 



(2) 



3 ^(L + \)T{N + L + 1) 1 — 4iV 2 



(63) 



4 r(n + l)r(L + l)r(2L + 2) (2L - 1)(2L + 3) 
These evaluations are carried out using standard representations and properties of the beta 
function and the hypergeometric function. 

Working back from these results to the original expression (f48|) for the energy correction, 
we finally have 



A4 = -a 5 ^,, j ° 



3vriV 3 
3 



' M 3 



2iV 5 (2L + 1) 



A 



1 1 - 4iV 2 

2 + (2L + 3)(2L- 1) 



(64) 



77T M 4 

The reader is reminded that this is the energy correction to order a e due to the correction 
AR = -R(p 2 ) — Re, and that its calculation used first-order perturbation theory only. 



D. Correction IRIRI + IRIRIRI + 



Finally we address the rest of the series for U, given by (p5|). 

It will be proven below that to the required order R(p 2 ) can be approximated by its 
lowest-order value c/2/i, where the dimensionless function of the masses c is defined by (p6|). 
Then the series (|35|) becomes to lowest order 



C/ 6 °(r) = I(c/2n)I(c/2fj,)I + 7(c/2//)/(c/2 At )/(c/2 A t)/ + 
a 3 (c/2/i) 2 



(65) 
(66) 



r 2 (r + ac/2/i) 

The energy correction up to order a 6 is given by the expectation value of (|66|) over non- 
relativistic Coulomb wavefunctions. 

For L > 1 the lowest-order expectation value is just dictated by the expectation value of 
1/r 3 . It is 



AE, 



L>1 



(U G (r)) L>1 = -aV: 



(67) 



2N 3 L(L + l)(2L + 1) 

For L = a logarithm occurs. Unfortunately expectation values of quantities of the 



form (|66|) over Coulomb wavefunctions do not seem to be listed in standard references. By 
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expanding the polynomial in the radial wavefunctions the result can easily be found as a 
double sum from to n. But in case this kind of integral arises later in real atomic two-body 
systems, at the cost of a few more lines of calculation we present an evaluation in closed 
form. 

In terms of the variable z = 2f3 Q r = 2a fir /N, the radial wavefunction is known to 
be z L e- z / 2 F(-n,2L + 2, z) up to normalisation, where F is the confluent hypergeometric 
function. Recalling that L = 0, and defining a = a 2 c/N, we need to evaluate 

r°° 1 o 
K a = e - z ^[F(-n,2,z)] 2 dz (68) 
Jo z + a 



While this expression has not been found in any reference, Landau and Lifshitz [jnj give 
an evaluation of a related integral: 

r°° 

J v = / e-'z u - 1 [F(-n,'y,z)] 2 dz (69) 
Jo 

T{u)n\ " (-n)^-7 + l),(7-^ (?(]) 
(7)n to I^M 1 ] 

where as usual = k(k + 1) • • • (n + k — 1), (k) = 1. Here the quantity 7 will be 2 since 

L = 0. The value v = 3 gives the normalisation integral J3 = 2. 

In the parameter a is small and a log(l/a) term will dominate, followed by 0(1) 
terms. These are the terms needed to find the energy to order a 6 . Also, if v is taken to be 
small in J v a 1/V term will dominate and the next terms will be 0(1). If [F{— n, 2,z)] 2 is 
written as 1 + B(z), in B(z) the lowest power of z will be one. Then it is clear that in both 
K a and J u , the leading contribution of the part B(z) is of order 1. Furthermore, those 0(1) 
terms can be calculated by replacing l/(z + a) by 1/z and z u ~ l by 1/z respectively. In other 
words, the contribution of B(z) to the 0(1) terms of both integrals is identical. In addition, 
the contribution of the "1" in each integral can be calculated explcitly. 

Therefore the way to evaluate the unknown integral (|68| ) up to 0(1) for small a is the 
following, (i) From the result ( f70| ) subtract the explicitly calculated contribution to J u of 
the "1" in [F] 2 . (ii) Evaluate the now leading 0(1) part of the result, neglecting terms of 
order v. (iii) Explicitly calculate the contribution of the "1" to K a to orders log(l/a) and 
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1, neglecting terms of order a. (iv) Add that to the result (ii), finally obtaining the log(l/a) 
and 0(1) terms of K a as required. The result is 

K a = log l - - (n + 2)C + - (n + l)V(n + 1) + 0(a) (71) 

and hence to order a 6 



L=0 



-V^flog^-^+DC+^-^w} ,72) 

in which iV = n + 1, V ; (- 2 ) = rflog[(r(2)]/dz and C is Euler's constant 0.5772 . . .. 

At the beginning of this Subsection we stated that it was sufficiently accurate to replace 
R{p 2 ) by cj2\x in the series (|3~5"D. Now it is necessary to prove that the correction to this 
approximation is of order a 7 or greater. With the definition 5R = R{p 2 ) — c/2/x, the 
correction to the a 6 term IRIRI will be 

I6RI(c/2fi)I + I(c/2fi)I6RI (73) 

We must prove that this term is of order a 7 only. Since Re — c/2/z = 0(a 2 ), we have 
5R = AR + 0(ct 2 ), and it is sufficient to prove that 

IARI(c/2fjt)I + I(c/2fj)I ARI (74) 

is of order a 7 . 

The calculations of Subsection C suggest that ( [71D should be of order a 7 , since (I ARI) 
was of order a 5 while (J(c/2/i)7) is of order a 4 . However, it is not quite obvious, since 
Ai2(oo) = —Re ~ — c/2/x, and therefore the rules of power-counting would permit ( |74|) to 
be of order a 6 . To see whether it is or not needs an explicit calculation. 

In momentum space, each term of (71) contains a factor of the form 



d\ e 2 e 2 

rA^l 2 )-^-, (75) 



(2tt)3( p -1)2 v '(l_q)a 

(preceded or succeeded by other factors) while in the first (a 6 ) term of the corresponding 
factor is 
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e 2 c e 2 



(76) 



(27r)3(p-l)2 2/U (l-q) 2 
We will exhibit counterpart terms in these two expressions and show that the one in (173) is 
smaller than the other by a factor 0(a). 

Again it is necessary to use the Schwinger unit sphere. With the mapping p — ► £ as 
before, as well as 1 — > rj and q — ► £, it is true that 



1-6 



1-774 



(P " 1) S 



(e-^) 2 
N' 



(77) 
(78) 



NLM 

Using the various relations given in Subsection C, expression (|76|) becomes for fixed L,M 



l-£ 4 v (2?r) 3 a 2 r . u (c/2 M ) y * m 1_C4 

2^ _/V' AT" JV ' iM ^ ' A L,N'N" Y N"LM IU " 







N'N' 1 



P 



(79) 



in which 



^KjV" = / ^F^, LM (77)^^—^^/(77) 



2/L* 1 — 774 

With the notations N' = n' + L + 1, N" = n" + L + 1, this is 



(80) 



.4 



(c/2 M ) 
L,N'N" 



2 L+1 r(L + i) 



n'!JV' 



27rr(Ar' + L + 1) V 27rr(iV" + L + 1) 



5 



(c/2 M ) 
L,N',N" 



[81) 



where, with y = 774, 



5 



/3c (i-i/ 2 )^g^(y)gy(y) 

z v v " " 2/, ./ , 1 ■ // 



(c/2 M ) 



(82) 



Expression (p2|) can be evaluated with the help of equation (|6l|). Supposing for definite- 
ness that n' > n", we easily find 

(c/ » _ /? Cv ^r(L + i)r(iV" + L + i) 



(83) 



L ' N '' N " 2/i n"!r(2L + 2)r(L + 1) 
Next we do the corresponding decomposition of (|75|) . The steps are the same, with AR 



replacing c/2/7. Using the representation (|53|) as before, we end up with 



m 3 M 3 



f3JL,n',n" 



(84) 



where 

•W = { V /I +^( 1 - vT^M^fo)} (85) 

The expression for Ji in ',n" is evaluated the same way that fl55|) was. Similarly to Ji in , it is 
found to be of order a 2 . 

Therefore, recalling that (3= 0(a), we finally have 

b£SIn» = 0(« 2 ) + 0(« 3 ) (86) 

while 

B%$%„ = 0(a) (87) 

Thus the expression (|75|) is of order a smaller than (f76[) , and so the correction term ([74) is a 
factor a smaller than the original quantity I(c/2fj,)I(c/2fi)I which is of order a 6 . It follows 
that the sum (|66"D does correctly represent the series (|35|) up to and including order a 6 . 

E. Summary 

The complete energy corrections of order a 5 , a 6 log(l/a), and a 6 , to the Bethe-Salpeter 
equation @ with a scalar Coulomb kernel, are given by the sum of: 

• Equation (37), from the two-body Bohr-Sommerfeld formula ([32]); 

• Equation (p3), from the correction IARI to U; 



Equation ( |67|) for L > 1, and equation (^|) for L = 0, which came from the sum (|66] 



of the rest of the Coulomb series for U. 



These corrections were calculated for every bound state, not just for a few low-lying 
states. The first correction was calculated by simple algebra. The latter two were calculated 
from first-order perturbation theory. As far as we know, calculations like these are beyond 
the power of the Salpeter reduction. 
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VI. CONCLUSION 



Advantages of the formalism shown here are that the local wave equation is derived 
directly from the Bethe-Salpeter equation, and that the two-body wave equation is accurate 
enough to predict the energy levels of the Bethe-Salpeter equation to first relativistic order 
correctly 

Two-body relativistic wave equations such as (B^) are easier to solve than Salpeter 
equations such as (^), because they are local in configuration space. The non-local 
components of the Salpeter equation, such as the free-particle kinetic-energy operator 
y/p' 2 + m 2 + Vp 2 + M 2 , and the correction factor Z p shown in equation (0), are replaced in 
the two-body relativistic wave equation by constants, such as E = \Jm 2 — f3 2 + \JM 2 — j3 2 , 
Z E — (rn+ M)/E, and the correction factor R E shown in equation ([19]). 

In addition, when the kernel of the scalar Bethe-Salpeter equation is a scalar Coulomb 
potential, the relativistic two-body bound-state wave equation is soluble exactly. The resul- 
tant two-body Bohr-Sommerfeld formula not only predicts the Bethe-Salpeter energy levels 
to order a 4 correctly, but it also allows much of the 0(a 6 ) correction to be evaluated by sim- 
ple algebra. Furthermore, the remaining a 5 and a 6 corrections are calculated by first-order 
perturbation theory only. No second-order perturbation calculations are needed. 

The results shown above suggest that it may be easier to solve Bethe-Salpeter equations, 
at least those whose binding interaction is a Coulomb potential, with the aid of a local wave 
equation rather than a Salpeter equation. It is hoped that the formalism developed here 
can be adapted to real two-body atoms such as the hydrogen atom and positronium. 

I thank the trustees of Springfield Technical Community College for a half-time sabbatical 
leave which assisted the completion of this work. I am grateful to G. Adkins for extensive 
comments. 
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